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Abstract 

The boundary effects on the Bose-Einstein condensation with a 
nonvanishing chemical potential on an ultra-static space-time are stud¬ 
ied. High temperature regime, which is the relevant regime for the 
relativistic gas, is studied through the heat kernel expansion for both 
Dirichlet and Neumann boundary conditions. The high temperature 
expansion in the presence of a chemical potential is generated via 
the Mellin transform method as applied to the harmonic sums repre¬ 
senting the free energy and the depletion coefficient. The effects of 
boundary conditions on the relation between the depletion coefficient 
and the temperature are analyzed. Both charged and neutral bosons 
are considered. 
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1 Introduction 


In this paper we will derive finite size effects on the Bose-Einstein condensa¬ 
tion of a Bose gas on an ultra-static space time R x M in the limit /3m — )■ 0. 
The Bose gas will be confined in a box B of arbitrary shape with volume V 
and surface area A, and it will be subject to either Dirichlet or Neumann 
boundary conditions. We will consider both charged and neutral gases. The 
metric on R x M is 

ds^ = —dt^ -I- 'jijdx^dx^. (1) 

Clearly the Minkowski space is a special case. 

The first step in our analysis will be to express the free energy A and 
the depletion coefficient W for the condensate in the presence of a non¬ 
vanishing chemical potential /i 7^ 0 as harmonic sums i.e. as sums of the 
form f{kl3m), involving the trace of the heat kernel. It is well known 
that in the relativistic regime it is the small values of /3m which is relevant 
for the Bose-Einstein condensation. The natural way to develop the asymp¬ 
totic expansion of a harmonic sum is to take the Mellin transform, determine 
the poles and residues of an appropriate meromorphic extension of the trans¬ 
form, and finally take the inverse transform (Mellin-Barnes integral) of the 
singular terms [1-3]. The next step will be the derivation of the small /3m 
asymptotic expansion (high temperature expansion) of A and Nf, by applying 
the Mellin transform techniques to the harmonic sums representing A and 
Nf,. No use will made of zeta function regularization which is the basis of 
the conventional approach to the high temperature expansion [4-11]. Thus 
our analysis in the present work will provide an alternative way of regu¬ 
larizing and developing the high temperature expansion. Next, we will use 
this high temperature expansion to study the effects of boundaries on the 
Bose-Einstein condensation on a manifold. It will be shown that in three 
dimensions and in the charged case the leading finite size correction is log¬ 
arithmic in /3m; whereas in the neutral case it will be a power law. Higher 
dimensional case will also be discussed. Our results for the free energy will be 
in accord with the ones derived in [8,10,11]. The usual objection against the 
use of Mellin transform in the case /x 7^ 0 (see e.g. Sec. 5.4 of [12] ) is a po¬ 
tential divergence caused by exponentially increasing factor of appearing 
in the transform. However we will show that large time decay properties of 
the heat kernel [13] imply that in fact the Mellin transform remains finite and 
the high temperature expansion proceeds without any divergence problem. 

Special relativistic Bose-Einstein condensation in thermodynamic limit 
was discussed by many authors [14-25], the interacting case in fiat space 
was studied in [26], and the curved space generalization was given in [27]. In 
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curved space, the relation between condensation and symmetry breaking was 
discussed in [10,11,37,38]. The hnite size effects in the non-relativistic flat 
case were studied in [28-34], those in the special relativistic case in [35,36,40], 
and those in the non-relativistic curved case in [39]. Finite size effects in 
relativistic gas of neutral bosons were presented in [41,42]. A discussion 
of the Bose gas inside a harmonic trap through Mellin transform can be 
found in [43]. For an approach to the Bose-Einstein condensation for both 
noninteracting and interacting systems on a manifold based on the global 
heat kernels estimates of [44] see [45]. 

Here is the brief summary of the present paper. In Sec. 2 we shall dis¬ 
cuss the relation between the asymptotic expansions of harmonic sums and 
Mellin transforms and shall work out the high temperature expansion of the 
free energy for non-vanishing chemical potential using this method. In Sec. 
3 we shall analyze the net charge of excited bosons using similar techniques 
and derive the hnite size effects on the relation between temperature and the 
charge density. In Sec. 4 we shall consider the free energy of the neutral 
Bose gas along the same lines, and in Sec 5. we shall investigate the effects 
of boundaries on the relation between temperature and number density of 
excited neutral bosons. In Appendix A we shall derive the asymptotic expan¬ 
sion of a certain function that is of fundamental importance in our analysis. 
In Appendix B we shall discuss further technicalities involved in the use of 
the Mellin transform to generate the high temperature asymptotic expansion. 

2 Free Energy of Charged Bosons 

The Klein-Gordon equation (KG) on R x M with external potential V is 

{—D + m'^ + U)^jJ = 0, ( 2 ) 

where 

□ = (3) 

Separating the time derivative we can write KG as 

- ^ ( 4 ) 

where 

A = (5) 

V' 

is the Laplacian on M. We will dehne the single particle Hamiltonian as 

H = -A + U + m\ (6) 
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We will assume that —A + f/ is a non-negative operator. This clearly implies 
H is positive. If the former operator fails to be non-negative we will still 
assume that H is positive. In the latter case without loss of generality we 
can redehne U and m? in such a way that the new —A -|- is non-negative 
and the new w? is positive. 

Free Energy of the charged boson is; 

A = — l^log j - 1 - log j 

^ a 

(^) 

Here A^r’s are the eigenvalues of —A -|- U which by assumption are non¬ 
negative. 

After using 

OO ^ 

- Iog(l -x) = J2^ (8) 

k=l 

and the subordination identity 






dt 


e 


A 

4t e 


—tx 


(9) 


in (7), and changing the integration variable to t = we get the free 

energy of the system as; 


A = ~ —= cosh(/c/3p) Fikjdm) (10) 

VTT 
k=l '' 

where 

F{x) = ^e^e-*(Tre^(^-^)). (11) 

The right hand side of (10) is a series of the form 

CXD 

=-^^g{kx), ( 12 ) 

where 

g(x) = coshi^x) F(x), x = fdm, ^ =—■ (13) 

m 

Such series are known as harmonic sums and their x —?■ 0 asymptotics can be 
determined from their Mellin transforms. The Mellin transform of a function 
f{x) is dehned as 


f{s) = {Mf){s) 


dxx^ f{x). 


(14) 
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A straightforward calculation shows that 


777 

(>lA)(s) = - C{s){Mg){s). (15) 

VTT 

Here C(s) is the Riemann zeta function which is a meromorphic function on 
the whole complex plane with a single simple pole at s = 1. 

In general the Mellin transform (Al/)(s) of a function / is a holomorphic 
function of s in a vertical strip (a, b) bounded by the vertical lines s = a 
and s = b. Largest such strip is called the fundamental strip of the Mellin 
transform. The small x asymptotic of f{x) is determined by the poles and 
residues of a meromorphic extension (if it exists) of the Mellin transform 
beyond the fundamental strip. This is done via the inverse Mellin (Mellin- 
Barnes) transform, 


{M-^f){x) 


rc+ioo 


ds f{s)x^ 


( 16 ) 


where the vertical line s = c lies in the fundamental strip. One has to 
complete this contour into a closed rectangular one enclosing the poles of 
the meromorphic extension of {Aif){s), then performing the integral one 
gets the asymptotic expansion of f{x). The details of the method can be 
found in [1-3]. The correspondence between the poles of {Aif){s) and the 
asymptotics of f{x) is as follows. If the Mellin transform has the singular 
expansion with residues A{w, k) 


w,k 


A(w, k) 
(s — 


where x refers to the singular part of a function, then 

w,k 


(17) 


( 18 ) 


Notice the presence of logarithmic terms which arise from multiple poles. 
In order for this to be really an asymptotic expansion one must make sure 
that the contributions coming from the extra horizontal and vertical lines 
used to close the contour are negligible. There are several different suffi¬ 
cient conditions which guarantee this. A sufficient condition suitable for our 
purposes will be 

/ OO ^ 

dr \f{a + iT)\ < oo. (19) 

■OD 
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Here a is a real number lying between two consecutive poles of f{s) with dis¬ 
tinct real parts. This part of the problem will be discussed in the Appendix. 

Before we apply this general method to (64) we must determine the x —?■ 0 
singular asymptotic of F{x), as will be shown below this will be crucial in 
the determination of the poles of (A4A)(s). The derivation of this singular 
asymptotics will be given in the Appendix A, here we just quote the result: 

4 

F{x) X ^ + b' log(x/2). (20) 

n=l 

Here the coefficients b are given in terms of the heat kernel coefficients of 
If 



^ OO 

k=0 

(21) 

&-4 

= 16m^ao 

(22) 

^-3 

= Ay/%rri^ai/2 

(23) 

b-2 

= 4(—-l- mai) 

(24) 

b-1 

= 2[-\pKvr?axi2F \px 

(25) 

b' 

= — m^oo -l- 2mai — m~^a2- 

(26) 


The behavior of F{x) for large x will also be needed in the following dis¬ 
cussion. Let A denote the hrst nonzero eigenvalue of —A subject to either 
Dirichlet of Neumann boundary conditions. If the smallest eigenvalue of H 
is nonzero then according to the Lemma 14.1 of [13] we have 

Tre-^* < Tr(e--^*°) for t > to- (27) 


Using the heat kernel asymptotics 


Tre-^^° 


U 


(47rto)^^^ ’ 

we see that by choosing to small enough we get 

3 V 


io —t 0 


(28) 


Tre 


-Ht 


< 


2 (47rot)^A 


(29) 
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Using this bound in (11) we obtain the bound 


F{x) < 


V 3^2 




X = 


V 3^/2 





(47rto)^F X (47rtQ)3/2 x 2x 

Here 7^1/2 (x) is the modihed Bessel function of order 1/2. Here we used 


(30) 


Ki/2{x) = 



(31) 


Thus 

F{x) = O [x~^^'^e~^) , X —)■ cxo. (32) 

If the smallest eigenvalue of H is vanishing, as in the case of Neumann 
problem for the Laplacian, we again have from Lemma 14.1 of [13]: 

Tre-^* < 1 + (Tre-^*° - l) (33) 


Proceeding as above we again get (30) and (32). 

We will see shortly that (Mg)(s) is holomorphic for 3fJs > 4. Now the 
main problem is the presence of exponentially increasing cosh term which 
makes it difficult to meromorphically continue (j\4g)(s) to 3fJs < 4. In order 
to hud the meromorphic extension we will write A as 

OO 

.4 = -^We-“*"9„(fa) (34) 

where 

ga(x) = cosh(^x)F(x). (35) 

Here a is an arbitrary parameter with 0 < a < 1 — which will be used 
as a regularizing parameter in the integrals arising in the derivation of the 
asymptotic expansion. 

Using (20) we see that as x —?■ 0 


9a(x) ~ (i_4X + c_3X ^ + d_2X ^ + (i_iX ^ + do + • • • (36) 

and using (32) we see that as x —)• oo 

g^(x) = o (x-3/2 g-(i-s-a)x^ ^ ^ 37 ) 
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where the coefficients d are related to the coefficients b. The hrst few are 


d_4 = 6_4 

d-3 = b-3 + a&_4 

d-2 = b-2 + a6_3 + 




b-A 


j 1.1. + (a^ + 3a^^) 

d—\ — 6—1 + otb—2 H-;;- 6—3 H-;;- 6—4 


do — cr6_i T 


2 6 


-6-3 + 
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Let’s also recall the hrst few heat kernel coefficients [47] 

1 / 

Oo — 


fll/2 — 


Oi = 


(47r)3/2 

A 

r - 

IGtt 

1 

6(47r)^/2 


dfig{—6U — R) + 


2K 


'B 


IdB 


(38) 

(39) 

(40) 

(41) 

6-4 

(42) 

(43) 

(44) 

(45) 


Here R is the Rieniann cnrvatnre of M, K is the extrinsic cnrvatnre of the 
snrface dB of the box, and dfig is the Rieniann volnnie element on M. In 
ai/ 2 , T = 1, — 1 for Nenmann and Dirichlet bonndary conditions, respectively. 
Now the Mellin transform of A is given by 


-m 


{MA){s) = -^C{s){CMga){s,a). 
Here CAi is the Laplace-Mellin transform [46] 

poo 

(£j\4ga)(s,a) = / dx ga{x). 


(46) 


(47) 


From the behavior of the integrand for large x (37) we conclnde that the 
npper limit of integration is not problematic. On the other hand for a > 0, 
g{x) = 0{x~'^) as x —)■ 0 and we see that the integral is convergent for > 4. 
Write {CAiga){s,a) as 


{CAiga){s,a) = / dxx^ \ga{x) — d_4X "^] — 


dxx^ 


'’d_4X 


dxx^ [ga{x) — d-4X '^] + d_ 4 a'‘ ®r(s —4). (48) 














From (36) we get ga{,x) — = 0{x~^) and we see that the lower limit 

of the integral is now hnite for > 3. On the other hand from (37) we see 
that the upper limit of the integral is also hnite. From the second integral 
in the above expression we get a pole at s = 4. Thus we get a meromorphic 
continuation of {CAig){s, a) to the region 3f?s > 3. Continuing in this manner 
and subtracting more and more terms of the asymptotic series (36) from g{x) 
we get a meromorphic extension of {CAig){s,a) into the half-plane to the 
left of = 4. The resulting expansion is 

{CMga){s, a) ~ d_ 4 a^-T(s - 4) + d_ 3 a^-T(s - 3) + d- 2 a^-T(s - 2) 

+d_ia^-T(s - 1) + doa-^r(s)... (49) 


Note that all the poles are simple. The pole at s = 4 comes from the d-^ 
term, while the one at s = 4 comes from both (i _4 and d^ terms, etc. 

Coming back to (A4A)(s) = —rmT~^^‘^C{s){CAig){s, a) we see that the 
poles are the same as above except that the pole at s = 1 is now a double 
pole because of the Riemann (^-function which behaves as 

C(s) ^-^ + 7 , (50) 

s — 1 

Here 7 is the Euler-Mascheroni constant. On the other hand we also have 
the singular expansion 


F(s 


(- 1 )"* 1 
ml s — n + m 

m=0 



(51) 


Using these we get the behaviour of (A4H)(s) around its poles and the cor¬ 
responding asymptotic expansion of A. 


Let us write 


A ~ -F ... 

(52) 

where is the contribution from the pole at s = m. 


For s = 4 


(A4H)(s) X ^C(4)5-4-T 

s — 4 

(53) 

Thus the corresponding asymptotic is 


1 

1 

-0 

07 

(54) 

For s = 3 

...... . -m C(3) , 

{MA){s)-^^b_s 

VTT S — 4 

(55) 
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Thus the corresponding asymptotic is 


dl(3) = ^C(3)^-3(/3"i) 

/n 


-3 


For s = 2 

( KA 4 V 

-^/F s — 2 2 

Thus the corresponding asymptotic is 






-2 


(56) 


(57) 


(58) 


For s = 1 


(-M^)(5) 


—m ( 1 


TT V s — 1 


+ 7) I -^-P('^-^5')(s)|s^.i + 


^6_3 + 6_i 


s — 1 


m 

¥ 


V ^-3 + ^-1 


(5-1)^ 


FP{CMg){s = 1,a) + 7 [f 6.3 + 6 i 


s — 1 


Here FP means the hnite part of the expression on the right 


(59) 


FP(CAig){s = 1, a) = {CM.g){s = l,a) 


^6_3 + b_i 


s — 1 


dxx"^ ^ cosh(,^a;)F(a;)— 


^&-3 + &-1 


s — 1 


(60) 


So 


^(1) = JIL 

'TT 


>2 

yfc-3 + &-1 


—m 

H—^ I 7 


y fe-3 + fe-i 


(/3m) ^ log(/3m) ^ 

+ FP(£>l^)(l,a) ) (/3m)-\ 


( 61 ) 
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Finally expressing all the b coefficients in terms of heat kernel coefficients 
throngh (38) and (22) we get: 


= 


-16 


TT 


C(4)m^ao(/Sm 


,-4 


^(3) = —4C(3)m^ai/2(/3r?7,) 

4m^ 


-3 


= 


c(2) [(2yU^ - m^)ao + oi] {/3m] 


-2 


= [ 2 (/i^m - m^)ai/2 + 2ma3/2)] (/ 9 m) Mog(/ 9 m) ^ 

( Tfl \ 

7 [2(p^m - m^)ai /2 + 2ma3/2] + -^FP{CMg){l,a) j (/9m)“^ 

(62) 

Note that the resnlt does not depend on the regnlarization parameter a. This 
resnlt is in accord with [8,10,11]. 

The above analysis can be repeated in dimension d > 3 with the result 


A 


—((d + 1) 


(2m 


id+l 


r 


TT 


d + 1 


ao {(3m) 


— (d+l) 


- C{d) 


(2m) 


r (^0 ai /2 {13m) 


—d 


c(d - 1) r 


TT 


d- 1 


([—m^ + (d — 1) /i^] ao + Oi) {(3m) 

(63) 


3 Charge 


Using the subordination identity (9) the net charge of the excited particles 
can be written as 


Qe 


E 

(7^0 


1 


1 


g^(Ao--l-m2)i/2-^^ _ g/3(Ao-+m2)i/2+/3^ _ 




k=l 


(64) 

(65) 
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Here 

rli 2 

^e^e-*(Tr'e^^), (66) 

where Tr' means the ground state is omitted in the trace 

Tr' = Tr e*^ - (67) 

Here Aq is the lowest eigenvalue of —A. Note that the heat kernel is now 
given as 

CX) 

Tr'e*"^ ~ (68) 

k=0 

where a 2 i_^_^ = a 2 i +3 — (—l)^Ao//!, and a'^ = a^, otherwise. The small x 
asymptotic expansion of Fi is still given by (22) with the only difference that 
instead of {a^} coefficients we now have {a),} coefficients. Moreover, since 
Fi{x) < F{x) the estimates (30) and (32) for F{x) are also valid for Fi{x). 
The high temperature expansion can now be obtained by starting with 


the expression 



(69) 

where 

fa{x) = e"'^a;sinh(,^x)F(x), 

(70) 

and proceeding as in the case of free energy. The small x asymptotic of fa 
is obtained by combining the asymptotic of F{x) with the Taylor expansions 
of sinh(^x) and e"*: 

faix) ~ C- 2 X ^ + C-iX ^ + Co + . . . 

(71) 

where the hrst few c coefficients are given in terms of b coefficients as 

C-2 = ^b-4 

(72) 

c-i = + a^b _4 

(73) 

Co — ^b_2 + oi^b_3 2 6 ) 

(74) 

On the other hand using (32), for a; —)■ cxo we get 



(75) 

Now the Mellin transform of Qe is given by 


{MQ,){s) = ^as){CMfa){s,a). 

Vvr 

(76) 
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roo roo 

{CAifa){s,a) = / dx[f (x) — c_ 2 X~‘^] + / dxx^~^ e~°‘^c_ 2 X~‘^ 

Jo Jo 

dxx^~^ e~°^ [f{x) — C- 2 X~‘^] + C- 2 a‘^~^T{s — 2) 

(77) 

Proceeding as in the calculation of the free energy we get the expansion 
{CA4fa){s, a) ~ c_ 2 a!^“^r(s — 2) + c_iQ!^“^r(s — 1) + Coa“®r(s) + ... (78) 



Note that all the poles are simple. The pole at s = 2 comes from the c _2 
term, while the one at s = 1 comes from both c _2 and c = — 1 terms, etc. 
For s = 2: 

(MQ,)M X ic(2)Ki.-4) -2— (79) 

^/7^ s — 2 

Corresponding asymptotic is 

Q(2)=vr-'/2C(2)(eM(«"' (80) 


For s = 1: 

(MQ,)(s) ~ IT-'/" + 7 ) (^FP(£X/)(ii) 


TT 


1/2 


(5-1)2 


[7e6_3 + FP(£M/)(5 


l,a)] 


1 

s — 1 


(81) 


Finally for s = 0: 


(MQe)(s) X —C(0) 


TT 


— &-4 + ib-2 
0 


Thus we arrive at the high temperature expansion 

Qe ~ QT’ + +... 

with 

QP = 7r"^^\(2)^6-4(/3m)"^ 


(82) 


(83) 

(84) 


QP = -TT (^ 6 - 3 ) (/5m) Mog(/5m) ^ + 

+ 71 - 1/2 [^^^b.s) + FPiCMfa){l,a)] {/3m)-\ (85) 
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Using (22) we can express the above in terms of the heat kernel coefficients 
as 
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-2 


TT 


= - (4myuai/2) (/9m) ^ log(/9m) ^ + 


TT 


(/9m) 


-1 


In higher dimensions, i. e. for d > 3, a similar analysis gives 




TT 


d + 1 


oo (/9m 




C(d - 2 ) 4,, r (^1^ 


C(4-3)2!!^44,r 


TT 


d- 1 


( 86 ) 


2 (d — 1) 2 

—m H--—/i I oq 9“ 


(87) 


To investigate the effect of the bonndary on the critical temperatnre we 
start by determining T as a fnnction of iVo, the number of particles in the 
ground state, for a hxed number of particles N. Let us hrst consider the case 
d = 3. Taking into account the leading boundary contribution in (86) we 
arrive at the equation 


2U A 

-C(2)p.m^(/9m)“^ — b—^m{(3m)~^ log(/9m)“^ = N — Nq. (88) 


TT^ 


dvr 


In the case of condensation from 


Qo — 


1 


we get 


_ I g/3(m+/4) _ 
1 


/i = ±m — 


PqoV 


0{V 


-2^ 


(89) 


(90) 


Here and below go = Qo/V and q = Q/V. If we set y = ksT/m and divide 
both sides of (88) by V 


2 2 2 

—C{2)ymy - —/rm?/log?/= g - go- 
TT^ dvr 


( 91 ) 
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Here e = AjV is the small parameter in terms of which we can expand the 
solution perturbatively 

= 2/ = 2/0 + £2/1 + • • • • ( 92 ) 

m 

Now we see that the sign of g — go is controlled by the sign of /r; g > go for 
/X > 0 and g < go for /x < 0. The two cases can be worked out simultaneously. 
Since the error in replacing /x by ±m is 0{V~^) we have (here we are assuming 
that the geometry of our box B is regular in the sense that 1/V « A/V) 


'2, E 

—C(2)m^2/^ - r—m^y\ogy = |g - go 
TT^ AtT 

Then the solution is given as 


2/0 


2/1 


( 7r^(l9-9o|) \ 

V 2C(2)m3 ) 

1 ^"(l9-9ol) 
32mC(2) ^ 2C(2)m3 


(93) 


(94) 

(95) 


At this point it is tempting to set no = 0 to determine the critical tem¬ 
perature, however this is an inconsistent step for hnite V, since it implies 
a divergent /x, and is justihed only in the thermodynamic limit i.e. after 
taking the H —)■ cxo limit. This situation is similar to the non-relativistic case 
discussed in [39]. 

The analysis of the higher dimensional case is similar. Setting y = 
(/3m)e = A/V, /x = m and making the perturbative expansion (92) 
as in the three dimensional case, we get the results; 

1 

d-1 

(96) 


i\q - qo\) 

m2'^ C{d - l)T 


_ A C,{d-2) r(f) 
4m(d-1) C(d-1) r 

for d dimensional charged boson. 


(97) 


4 Free Energy of Neutral Bosons 

Starting, as in the charged case, with the subordination identity (9) we can 
write the free energy for a gas of neutral bosons 

A = ^Y^ [log (l - j (98) 
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as 


A = --^—e^X^Fik/Sm) 

CX) 

20F^ 


( 99 ) 


where 


ga{x) = ~ d-iX ^ + d-'iX ^ + ... x ^ 0. (100) 

In this case d coefficients in terms of b coefficients are; 

(i_4 = 6_4 

d -3 = 6-3 + 6-4(0! + 


d—2 — 6—2 + b—^ipi + ^) + 6—4 


(« + 0^ 


d-1 = 6-1 + 6-2(0 + 0 + 6-3 ^ + 6-4^^^ 

(o + £)2 (o + £)3 (o + £)^ 
d ^ = 6- i(o + 0 + 6-2 ^ , ^ +6-3 ^ +^-4 (101) 


6 
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Now the Mellin transform is 


—m 


where 


(A4/l)(s) = ^^((s)(£Mga)(s,a), 

2y 71 


(£A4g)(s,a)= / dxx^ e °‘^ga{x). 


( 102 ) 


(103) 


The singular asymptotic can be worked out exactly as in the charged case. 
Thus we have 


{£Mg){s, a) ~ d- 4 (o)^-T(s - 4) + d- 3 (o)^-T(s - 3) + .... (104) 

Now the poles are at s = 4, 3, 2,1,.... Again the pole at s = 4 comes from 
d-i term, the pole at s = 3 comes from both d-^ and d -3 terms, etc. We 
stress that the equation is in the same form as in the charged case. Let us 
express the free energy as; 

A ~ + A^^'> + + ... (105) 
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where 

= 


-m 


20F 


C(4)(16m^ao)(/Sm 


,-4 


^(3) = ZUL 
2^ 

A(2) = 


C(3)(16/im^ao + A^/^iw?al/2){|3m) 


-3 


= 


20F 


m 


20F 


C(2) [4(2m/i2 — rrV‘)aQ + A^/^iflmal /2 + 4mai] [fim] 


-2 


/i^ - 4 m^/i ) ao + 2 -v/ 7 r(/i^ - m^)ai/2 + 4 /iai + 20?03/2 


X (/9m) ^ log(/9m) 


\-i 


m / 
00 0 


-0 - 4m^/i ) ao + 2-v/7r(0 - m^)ai/2 + 4/iai + 20ra3/2 


+ FP(/:A4r?„)(l,a))(/3m)-^ 

Similarly, for higher dimensions, i.e. for d > 3, we get 

ao (/9m)-(‘^+^) 

'd 


(106) 




20F 


- cM^[2,r 


20F 


d + 1 


ao + r 


01/2 


(/3m 


, —d 


- C0-1) 


(2m) 


d-l 


20F 


r 


d -1 


[(-m^ + (d- 1)0) ao 


2/1 r 


d 


01/2 + Ol 


(/3m; 


-(d-i) 


(107) 
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5 Number of Neutral Bosons 

With the subordination identity (9) we get for a gas of neutral Bosons 


= E 


1 


cr^O 




f; 


k=l 


2\/i 


-| 

k=l 


where 


fot{.x) = ~ d_^x ^ + d_j,x ^ + ... x —)■ 0. 

Again the Mellin transform is 


( 108 ) 


(109) 


( 110 ) 


(lMiV,)(s) = ^({s)(CMf)is,a), 

ly TT 


( 111 ) 


where 

POO 

{CMfa}is,a)= / dx x^~^ f {x). (112) 

Jo 

The asymptotic can be worked out exactly as in the previous cases with the 
result 


iCMfa){s, -0 X d-4(a)^"T(s - 3) + d_3(a)2-T(s - 2) + .... (113) 

Now the poles are ats = 3,2,l,.... Again the pole at s = 3 comes from d _4 
term, the pole at s = 2 comes from both d _4 and d _3 terms, etc. The high 
temperature expansion is then given by 
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where 


JV(S) 

nP 

nw 


1 

1 

-1 

1 


C (3) (16m^ao) (/9m) 

C( 2 )( 16 /im^ao + 4 \/^m^ai/ 2 )(/ 9 m)“^ 

(4(2m/i^ — m^)ao + 4 \/^/U,mai /2 + 4mai) (/9m)“^ log(/9m)“^ 

[7 (4(2m/U^ — m^)ao + Ay/wiimai /2 + 4mai) + FP(£A4/a)(l, a)] (/9m)“^. 


(115) 


In higher dimensions, i.e. for d > 3, a. similar analysis yields 






TT 


+ C('9 — 1) 


+ C{d — 2) 


( 2 m 


.d-i r 


( 2 m 


TT 


i(i -2 


2 /ir 


aoiPm) 
d + 1 


oq + r 


d 


fll/2 


(/9m 


I — d+l 


TT 


d- 1 


[(-m^ + (d - l)/i^) ao+ 


2 /ur ^0 ai /2 +ai] {/3m] 


-d+2 


(116) 


As in the charged case we can determine the dependence of temperatnre 
on no (n = N/V and rie = Ne/V) using this high temperature expansion. 
Let us start with the case d = 3. We again set y = ksT/m, e = A/V and 
note that from 


No = 


1 


(117) 


we get 


1 


/i ~ m — 


l3noV' 


(118) 


Since the difference /i — m is of lower order than e = A/V we replace p by 
m and get 


^C(3)mV + C(2) 

TT"' 


1 3 I ^ 2 

—m ± —m 
7r2 Stt 


= n-no 


(119) 
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The perturbative expansion y = yo + eyi + ... yields 


and 


mvi 

3C{3)m^ 


71^ 


VoVi + 


C{2)y^o = 

n‘^{n - no) 


2C(2)m3 

C{2)m' 

^2 

Voyi ± 

Sti 


-Vo = 0. 


The last equation can easily be solved for yi in terms of yo 

.C(2)7r yo 


yi = T- 


8m 3C(3)|/o + 2C(2)' 


( 120 ) 


( 121 ) 


( 122 ) 


On the other hand the first equation being a cubic has an exact solution, 
which however is not very transparent. Instead we can solve it by the method 
of detailed balance [48], assuming small m. Let us write it as 

C(3)!/| + nmvl = ~"°f (123) 


Here rj is an expansion parameter which will be set equal to 1 at the end of the 
calculation. Now plugging in the perturbation expansion yo = Zo + rjzi + ... 
we get 


((C(3)) V2(n-no))^/^ 


C( 2 ) 

3C(3) •••• 


(124) 


By working few more terms of the perturbation series one can easily see that 
it is an expansion in powers of m and the neglected terms are in fact 0{m). 
Note that the leading term above is the result obtained in [41]. 

The analysis can be generalized to higher dimensions, with the result 


{n - no) Cjd-l) 

mdC{d)r{^) dCid) 


(125) 


^ ^ r (f) _ C{d-l)yo _ 

4m T {^) [dC{d)yo +{d-l)C{d-l)y 

where d > 3. 


(126) 
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Appendix A: Asymptotic Expansion of F{x) 

We will now derive the asymptotic expansion, more precisely only the singular 
part of the asymptotic expansion, of the function F{x) as a: —)■ 0 and thus 
determine the coefficients (22). 

Write 


poo 2 

4 



/o 

n=0 


r°° 2 ^ 

'fo n=0 

Now split the hrst integral as 


—2\nl2 


(127) 


rl fCO 2 

4 

/ + / 


Jo Jl 

n=0 


(128) 


The hrst integral is absolutely convergent as x —?■ 0 since because of the 
heat kernel expansion the term in square brackets is 0(t~^) and therefore 
the integrand is as f —)■ 0. On the other hand in the second integral 

term and the boundedness of Tr'e^* as f —?■ cxd imply convergence. 
Thus we see that the singular asymptotic of F{x) as x —)■ 0 is given by 


F{x) X '^an/ 2 K{x), 


where 


FJx) = 


n=0 




(129) 


(130) 


.-0 ' 2 . 

Here i^’s are the modihed Bessel functions. Using the series expansion 
thereof for non-integer u [49] 


KJx) = 


71 csc(z/ 7 r) 


E 

.k=0 


T{k — z/ + l)fc! V2/ 


X\^k-v 


[v ^ -v] 


,(131) 
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and for integer n 




(132) 


collecting the relevant terms we get 




(133) 


From this we get the b coefficients as given in (22). 

Appendix B: Asymptotic Expansion Through 
Mellin Transform 

Here we will discuss the validity of the Mellin transform to generate the 
high temperature expansion. We will consider only the free energy for the 
charged bosons. The other cases can be treated similarly. Consider the Mellin 
transform {Aif){s) of a function f{x). In general {Aif){s) is holomorphic 
on a strip a < < b. It is assumed that {Aif){s) has a meromorphic 

extension into 3?s < a with poles at a = wi > W 2 > ■ ■ ■■ Let A{wr, K) be the 
residue with multiplicity of the pole at Wr- 
The inverse Mellin transform is given by 



(134) 


where a < c < b. The asymptotic expansion of f{x) is obtained by closing 
the above contour into a rectangular one where the added vertical contour 
at 3?s = (To < c lies between two consecutive poles, say Wr and Wr+i (as¬ 
suming the real parts of these poles are distinct).The expansion is derived 
by completing the contour in the inverse Mellin transform into a rectangular 
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contour enclosing the poles of the meromorphic extension of {Aif){s). For 
this to be an asymptotic expansion the contributions of the edges used to 
close the contour must be negligible. The following sufficient condition is 
suitable for our purposes. Suppose a is a real number lying between two 
consecutive poles of {Aif){s) and let us suppose 


dr \{Aif){a + fr)| < oo. 


(135) 


Then what we calculate is an asymptotic expansion of the function f{x). To 
summarize, if 


lim \{A4f){a + iT)\=0 for (Tq < a < c, 

\t\^oo 


(136) 


and 


dr\(Mf)(cro + st)| < oo, 


(137) 


then we hnd 


r kr 


fix) ~ K)X "'"(log 


X 


i kj' 


2=1 k 

We will apply this to our case where f = A, the free energy. 


(138) 


-m 


iMA){s) = —^C,is)iCMgo)is,a). 

'7T 


(139) 


Now |C(o' + ir)| increases algebraically as |r| oo [50] (provided a is not 
close to 0, which is true in our case). If we can show that for large |r|, 
\{CAiga)is, a)\ = 0(P(|r|)e“'^l'"l) where P is a polynomial and q > 0, then 
the hypothesis of the theorem will be satished. 

The meromorphic extension of {CAdga) we are interested in is 


dx 


xS-l^-ax 


-1 


9a[.X^ 


-E 

j=-4 


CjX^ 


+ 


-1 

E 

j=-4 


CjOt 


-3-s 


r(s + j). 


(140) 


Here s = a + ir = |s|e®'^, a = 1/2. For later convenience we also dehne 
s = s/|s| = p + iq. From [49] 

|r(a + ir)| = 0(|r|'"-^e-^’"l"l) (141) 


we see that the last four terms have the the desired property of fast decay. 
So let us concentrate on the hrst term. 
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Rotate our contour (0, oo) by 0 < 7r/2. Along this new contour C the 
integral is equal to 

r 

/ dz ... = l dxx'^~^e~°'^^[g{sx) — ^ CjX^]. (142) 

Jc Jo 

For real x let us write (11) as 

= 4j^ r ^ e-^e-"“Tre^“. (143) 

X ' Jo 

Using this expression to continue F analytically to > 0 we get the fol¬ 
lowing analytic continuation of g 

g^{z) = e“" cosh{^z)F{z). (144) 

We know that for real x, ga{x) — ~ ^(1)- Assuming the heat ker¬ 

nel expansion can be continued to complex time and repeating the derivation 
of Appendix A with x replaced by sx we get 


-1 



aisx) - CjF = 0(1). 

j=-A 

(145) 

Thus 


/ dz ... 

Jc 

{apY 

(146) 

Now 

S 

1 1 

SI _ 1'^ 1 


(147) 

Here we noted ref) = 

kll0l 

. So we have the desired fast decay property. 


We must also show that the contribution to the result of the big circular 
arc Cji of radius R that connects the original contour (0, cxo) to C is negligibly 
small. Setting z = we get 


\ga{z)\ < e“^“"®cosh(^Rcos0)|F(^)|. (148) 

Now notice that 

\F{z)\ < R-^>^F{Rcose). (149) 

Thus we obtain the bound 

\ga{z)\ < cosh.{^Rcos6)F{Rcos6). (150) 
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Using (32) in this we get 


\g^{z)\ <ci?-i/'e“^“^"cosh(ei?cos0)(/2cos0)-='/2e-^“^® (151) 

Since —l + a + .^ <0 and —l + a — <0 the prodnct of exponential terms 

is less than 1. Thus we arrive at 

Ifi'a(^) I < c(-Rcos 6')“^^^. (152) 


Now along Cr 


r 

r<t 

-1 

/ dz... 


\9a{z)\ + 

Jcr 

Jo 

i=-4 


< / de 

Jo 

This shows that the integral along Cr vanishes as i? —?■ oo. 


-1 


c{R cos + 

i=-4 


^ 153 ) 
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